ON GENERATING RELATIVE AND ABSOLUTE INVARIANTS 
OF LINEAR DIFFERENTIAL EQUATIONS 
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Abstract. A general expression for a relative invariant of a linear ordinary 
differential equations is given in terms of the fundamental semi-invariant and 
an absolute invariant. This result is used to established a number of properties 
of relative invariants, and it is explicitly shown how to generate fundamental 
sets of relative and absolute invariants of all orders for the general linear equa- 
tion. Explicit constructions are made for the linear ODE of order five. The 
approach used for the explicit determination of invariants is based on an in- 
finitesimal method. 



1. Introduction 

The first invariants found during the early days of development of the theory of 
invariants of differential equations were all relative invariants [BUI [3], and one of the 
very first determination of absolute invariants for differential equations is perhaps 
due to Brioschi [3], who obtained them as a quotient of two relative invariants. It has 
subsequently been shown that in the case of linear ordinary differential equations, 
every fundamental absolute invariant can be expressed as a rational function. 

These basic properties of relative invariants have led to significant progress in the 
study of invariants of differential equations and their applications, and such studies 
have been largely influenced by the work of Halphen [H[S], and Forsyth [BJ. However 
the methods used by Halphen, Forsyth, and earlier researchers on the subject were 
very intuitive and most often ad hoc methods requesting tedious calculations just 
for finding a couple of invariants. 

Based on recent advances in Lie group techniques [3 HI [9l [10] , infinitesimal meth- 
ods haven been increasingly used for the investigation of invariants of differential 
equations [HI H21 El HH [16] . But although these infinitesimal methods provide 
a more systematic route for the treatment of invariants of differential equations 
and transformation groups, they have been applied to the determination of rela- 
tive invariants of differential equations only in some very rare cases |17i 1 1 2 j , and 
even in those cases only some very specific relative invariants were obtained in 
the usual way as absolute invariants corresponding to partial structure-preserving 
transformations, in which some of the variables in the equation are kept constant. It 
therefore appears that a number of interesting properties of these relative invariants 
have not yet been uncovered. 

Restricting our attention to linear ordinary differential equations, we determine 
in this paper some properties of relative invariants, considered as semi-invariants 
of the full structure-preserving transformations of these equations. In particular 
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we show that every absolute invariant can be expressed as a quotient of a relative 
invariant and the fundamental semi- invariant, and we derive a general expression 
for these relative invariants. We show how these relative invariants can be used to 
obtain invariants of all orders via invariant differentiation. Our approach for any 
explicit determination of invariants is based on the infinitesimal method recently 
proposed in |10j , and which contrary to the former well-known method of [8] , does 
not require the knowledge of the structure-preserving transformations, but rather 
provides it. 

2. Basic properties of relative invariants 

Let G be a Lie group of point transformations of the form 

x = ((>(z,w;t), y = ip(z,w;r), (2.1) 

where r denotes collectively some arbitrary parameters specifying the group element 
in G. Consider on the other hand a family T> of differential equations of the general 
form 

n(x,y (n) ;/0=O, ( 2 - 2 ) 

in which yr n \ represents the dependent variable y = y(x) and all its derivatives 
up to the order n, and p represents collectively some arbitrary functions of x, 
or arbitrary constants specifying the family element in T>. We say that G is the 
equivalence group of (|2.2[) if it is the largest group of transformations that maps 
elements of T> into itself. In this case, the transformed equation takes the same form 
n(z, wi n \; ff) = 0, in terms of the transformed parameter 9, and (|2.1j) is called the 
structure-preserving transformations of (|2.2j) . By a well known result of Lie [18] . 
the transformations (12. ip induces another group of transformation G c acting on the 
parameter p of the differential equation, and we shall be interested in the invariants 
of this group action and its prolongations, and which are commonly referred to as 
the invariants (or differential invariants) of (|2.2[) . As sets, G and G c are equal, 
except that they act on different spaces, and both are often referred to as the 
equivalence group of (|2.2jl . Thus both sets will often be denoted simply by G. 

It should be noted that in (|2.2p , x and y each denote collectively all independent 
and dependent variables, respectively, so that the equation also includes in partic- 
ular all partial differential equations. However, we shall be interested in this paper 
in a linear ordinary differential equations of the general form 

yin) + 0iy (n-D + a2y (n-2) + . . . + ^ = ^ (3 3) 

where a,j — aj(x) are arbitrary functions of the independent variable x. The 
structure-preserving transformations of (|2.3p can be written in the form 

X = £(«), y = v(z)w, (2-4) 

where £ and r\ are arbitrary functions. Under (|2.4p , the transformed equation 
of (|2.3p takes the form 

+ A^-V + A 2 w {n - 2) + ■■■+ A n w = 0, (2.5) 

where the Aj — Aj(z) are the new coefficients. Let a denote collectively the co- 
efficients dj = aj(x). A differential function F = F(a,ar r -)), where ar r ) represents 
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as usual the derivative of a up to a certain order r, is called a relative invariant 
of (23]) if 

F(t ■ (a, a (r ))) = w(t) • F(a, a (r) ), (2.6) 

for all a and r, and for all t G G. In (|2.6p , the weight function w = w(t) must 
be a character of the group G. When w is identically equal to one, the function F 
is called an absolute invariant of (12. 3j) . For simplicity, an expression of the form 
F(r • (a, a( r ))) like in (|2.6j) will often be represented by F r , for a given function F. 

Let's assign to an expression of the form d k aj/dx k the weight j + k. We say 
that a polynomial function F in the coefficients aj and their derivatives has weight 
m if each of its terms has constant weight m. By combining a result of Forsyth [5] 
according to which all absolute invariants of (|2.3p can be obtained as rational func- 
tions, and certain results obtained by Halphen [5], we readily obtain the following 
result. 

Theorem 1. Equation (|2.3p has a fundamental set of absolute invariants consisting 
of rational functions, in which every element F — S1/S2 is the quotient of two 
relative invariants Si and S2 of the same weight m, each of which satisfies a relation 
of the form 

S(T-(a,a (r) )) =£(z) m S(a,a ir) ), that is S T = £(z) m S, (2.7) 

for every r G G, where S denotes any of the invariants S\ and S2 ■ 

The invariant F of Theorem[T]is also said to be of weight m. In general, a relative 
invariant S satisfying S T = £ r S for some r G IR is said to be of index r. We shall 
make the result of this theorem more precise in the next section. 



3. The fundamental relative invariant 

For a general group of transformations G acting on a manifold M, a semi- 
invariant usually refers to a function F satisfying a relation similar to that specified 
by (|2.6[) . and which has the simple form 

F(g-p)=w(g)-F(p), (3.1) 

for all g G G and p G M such that g ■ p is defined, and where the function w is a 
character of G. If we let v be a generic element in the generating system for the 
Lie algebra of G and denote by X the corresponding infinitesimal generator of the 
group action, then (|3. 1(1 implies that 

X ■ F = —XF, where (3.2a) 
A = dw(e)(v), (3.2b) 

and where dw(e) denotes the differential of w at the identity element e of G. A 
function F satisfying (|3.2|) is often called a X-semi-invariant. By a nontrivial semi- 
invariant, we shall mean a semi-invariant which is not an absolute one, and which 
in particular is not a constant function. We now establish the following result that 
relates every semi-invariant to a fundamental set of absolute invariants. 

Theorem 2. A function F defined on M is a X-semi-invariant of G if and only if 
it is of the form 

F = F $, (3.3) 
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where Fq is an arbitrarily chosen nontrivial X-semi-invariant, and $ is an absolute 
invariant, that is $ = $(Ji, . . . ,J a ), where {Ii,...,I s } is a fundamental set of 
absolute invariants of G. 

Proof. Suppose that in a coordinates system {u\ , . . . , u q } = u of M the infinitesimal 
generator X of G defined as in (|3.2p has the form 



X = X 1 d Ul +--- + X q d. 



q '-'Uq J 



where Xj = Xj(u). Then by a well-known result [19] . the equivalent system of 
characteristic equations associated with (|3.2[) is given by the sequence of q equalities 

du\ du2 du q dF 

Xl X2 X q —XF 

in which the first q — 1 equalities are the determining equations for the absolute 
invariants. To find the general solution to this system, we can first find the integral 
resulting from a combination of the last fraction dF/(—XF) with any suitable ones 
from among the first q fractions. We may assume without loss of generality that 
such a combination is given by an equation of the form 

dF —Xduk 



F X h 



for some k, 1 < k < q. 



This equation clearly has separable variables, since F is considered as a new in- 
dependent variable added to the coordinates system {u\, . . . , u q } . Its integral can 
be written in the form Fv = C q +i, for some function v = v(u). If we denote by 
Ij(u) = Cj, for j — 1, . . . , q — 1, the other integrals corresponding to the first q — I 
equalities in (|3.4p , and in which the Cj are arbitrary constants, then the Ij are the 
absolute invariants of G and the general solution of (|3.4p takes the form 

$i(ii, . . . ,I q ,Fv) = 0, or equivalently Fv = <&(ii, . . . ,I g ), 

for some arbitrary functions $1 and To obtain the function Fq of the theorem 
we only need to take Fq = 1/z/, and it readily follows from the Leibnitz property 
of the derivation X and the definition of an absolute invariant that Fq is a A-semi- 
invariant. □ 

It should be noted that Eq. (|3.3|) also holds with the same Fo for all prolonga- 
tions of the group G. Since _Fo is nontrivial, it follows in particular that the set 
{Ii, . . . , I q , Fo} is functionally independent. In contrast to the case of absolute 
invariants, not every function of semi-invariants is again a semi- invariant, and al- 
though the set of all semi-invariants of G forms a group under functions multiplica- 
tion, the sum or difference of two semi- invariants is not in general a semi-invariant. 

For simplicity, we let I denote collectively all elements in a fundamental set of 
absolute invariants of Eq. (|2.3p . We shall also often denote by to an element of the 
form (a, a( s )), which can be viewed as an element in the sth-jet space determined 
by the independent variable x and the dependent variable a of Eq. (|2.3j) . 



Corollary 1. Suppose that Si is a relative invariant of Eq. (|2.3p of index k. Then 
every relative invariant S2 of order m of the same equation can be put into the form 

S 2 = S™ /k 1>, (3.5) 
for a certain function $ = $(/). In particular S™ / S% is an absolute invariant. 
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Proof. According to Theorem[2j to prove the first part of the corollary, we only need 
to show that is a A-semi-invariant for the same function A as S 2 , and by (|3.2[) . 

we only need to show that S\^ is a relative invariant of the same index m as S*2. 
But since (Si) T = £ k S\ by assumption, it readily follows that (S™^ k ) T — £ m S™ l ^ k . 
For the second part we notice that since S p has index rp for every relative invariant 
S of index r, we must have (5J n /5^) r = S[ n /S%. □ 

The corollary says that every relative invariant of (|2.3|) can be expressed as a 
product of an arbitrarily chosen, but fixed relative invariant and an arbitrary func- 
tion of /. We call such a fixed relative invariant the fundamental relative invariant 
of (|2.3p . and denote it and its index by Sq and a, respectively. We can now give 
without any need of calculations a simple formula relating every absolute invariant 
with the fundamental relative invariant. 

Corollary 2. Eq. (|2.3|) has a fundamental set of absolute invariants in which every 
element has the form F — Si/Sq 1 ^ for some relative invariant S\ of index m. 

Proof. By Theorem [T] we may write every fundamental absolute invariant of index 
m in the form R\/ R 2 , where R\ and R2 are relative invariant of index m. It follows 
from Corollary □ that R 1 = S™ /a <S> 1 and R 2 = S™ /a <$> 2 , for some functions $1 and 
$2- The corollary is thus proved by taking Si = 5™^° $i/$ 2 - □ 

Corollary 3. The weight and the index of every relative invariant coincide. 

Proof. Suppose that the relative invariant S has weight \x and index m. Then 
S/S™^ is an absolute invariant whose weight is the index of S™^° , which is m by 
Theorem [1] By the same theorem, this weight is the same as the weight of S, 
and this proves the corollary. □ 

Corollary 4. 

(a) Suppose that {I\, . . . ,I p } is a fundamental set of absolute invariants of a 
certain group of equivalence transformations G, and that Fq is a nontriv- 
ial A- semi-invariant of G. Then, {Fq, FqIi, . . . , FqI p } is a maximal set of 
functionally independent X- semi-invariants of G. 

(b) If Ij — Sj/S™ 3 ^ a are the fundamental absolute invariants of Eq. (|2.3p for 
j = 1, . . . ,p, where rnj is the index of Sj, then 

s = {s™ /a ,s™ /rni ,...,s; n / m ?} 

is a fundamental set of relative invariants of index m of Eq. (|2.3p . 

Proof. Since Fq is a nontrivial semi-invariant of G, the given set {Fq , FqIi , . . . , FoI p } 
clearly forms a functionally independent set of semi-invariants of G, by Theorem [2] 
By the same theorem, any other semi- invariant of G has the form S = F <S>, where 
$ = $(/!,..., I p ), and this readily proves part (a) of the corollary. 

On the other hand, every element in the given set S is clearly a relative invariant 
of index m. if we replace each Ij by 

the resulting set formed by the J™/™ 3 [ s a fundamental set of absolute invariants of 
the same index m of Eq. (|2.3[) . Thus by (|3.2[) the elements of S are A-semi-invariants 
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corresponding to the same function A. Therefore, part (b) of the corollary readily 
follows from part (a). □ 

It also follows from part (b) of this corollary that {So, Si,..., S p } is a functionally 
independent set of relative invariants, but whose elements do not have the same 
index. 

4. Application to the determination of invariants 

It is well-known, by a result of Lie, that all higher order differential invariants 
of a transformation group can be obtained from a generating system of lower order 
ones by means of invariant differentiation. This often reduces the problem of deter- 
mination of invariants to finding a generating system of invariants and the invariant 
differential operators. To begin with, suppose that we know two functionally inde- 
pendent absolute invariants Iq and I\ of equation (|2.3|) , and that these are rational 
functions of the form 

Io = R°JS k , h = SIJSZ (4.1) 
where 5*0 is the fundamental relative invariant of index a already introduced, while 
i?o and Si are relative invariants of indices k and m, respectively. We also assume 
that S*i is of order /iasa function of the independent variable x, and that all other 
relative invariants are of order at most /x in x. Then by means of the differential 
operator C,D X , where £ = 1/Iq, Iq = dlo/dx and D x — d/dx, we can generate a 
new absolute invariant /* = (^D x (Ii), whose order is fi + 1 in general. In terms of 
the relative invariants in (14.11). we have 



r = h Ro (mSiSp ~ rSoS[) 
1 Io Si (kR Q S -rS Rd)' 1 ' ' 

where h! = dh/dx for every function h = h(x). It follows from Corollary [5] that 
to find an absolute invariant of higher order ji + 1, we only need to find a relative 
invariant of order /j, + 1. We now introduce the notation 

^(Ri,R 2 ) = m 1 R 1 R' 2 -m 2 R 2 R' 1 , x^Jfa) = (4.3a) 

V0OR1.-R2) = Xo(Ri,R2) = RT 2 /RT (4.3b) 

for any relative invariants Ri and R2 of respective index mi and m,2- For simplicity 
of notation, when the second argument R 2 is fixed and there is no possibility of 
confusion, we shall set <p(Ri) — (p(Ri,R2) and x(-^i) = xi^ii^a)- By multiplying 
I* by the absolute invariant I0/I1 and the relative invariant Si/Rq, we obtain 
a relative invariant of the form F = tp(Si, So)/ip(Ro, So). By an application of 
Theorem^ it is easy to see that both ip(Si) and <f(Ro) are relative invariants, and 
they clearly have indices m + r + 1 and k + r + 1, respectively. Moreover, tp(Si) has 
the required order ji + 1, and it gives rise to the absolute invariant 

X (Si,So)= {m+r+iy (4-4) 

Owing to the arbitrariness of the function Si , Eqs. (|4.3[) and (|4.4[) show that starting 
with any relative invariant Si of order fx, we can construct an indefinite sequence 
<Pq(Si) — ip q (Si) of relative invariants and Xq(Si) = x q (Si) of absolute invariants, 
each having a general term of order fi + q. To write down the expression for 



ON GENERATING RELATIVE AND ABSOLUTE INVARIANTS 



7 



the general terms of these sequences, we only need to note that (p q (Si) has index 
6(q) = m + q(a + 1). Consequently, we have 

<p q (S 1 )=9(q-l)<p<- 1 (S 1 )S' -rS D <p 9 - 1 (S i y t Xq (S 1 )= [ ^^f , (4.5) 
where 

= d( f p*- 1 {Si))/dx. 

A similar sequence of relative and absolute invariants can be obtained using the 
relative invariant (f{Ro, So). 

In the case of ODEs only one differential operator of the form C,D X is required and 
once this operator is known, finding fundamental sets of invariants only requires a 
suitable generating system of invariants, and these invariants must be explicitly cal- 
culated. As already mentioned, earlier methods for finding these invariants such as 
those used in [5] were very intuitive and quite informal, and they were designed only 
for linear equations. However, based on ideas suggested by Lie [20], and starting 
with some works undertaken by Ovsyannikov [7j, the development of infinitesimal 
methods started to grow and a number of Lie groups techniques for invariant func- 
tions have been proposed in the recent scientific literature [H [2TJ \§\ [TO] • These 
Lie groups methods provide a systematic means for finding invariant functions and 
some of them have far reaching immediate applications (see [2TJ[10]). For instance, 
the method of [10] that we shall use for the explicit determination of the invariants 
also provides, in infinitesimal form, the structure-preserving transformations of any 
differential equation. When it is applied to Eq. (|2.3p , the coefficients a,j (x) of the 
equation are considered as dependent variables on the same footing as y and the 
infinitesimal generators X° of the equivalence group G c takes the form 

n 

*° = /fl» + X>0««. (4.6) 
»=i 

The explicit expression of X° depends on various canonical forms adopted for 
the general linear ODE. By a change of variable of the form x = z and y — 
exp (— J a\dz\ w, Eq. (|2.3p can be put in the form 

iflW + b 2 w (n ~ 2) +■■■ + b n w = (4.7) 

which is deprived of the term of second highest order, and in which the bj = bj(z) 
are the new coefficients. Similarly, by a change of variables of the form 

{z,x} = — — TTT a 2, V = ex P ( - / a-idz ) w, (4.8a) 

n(n-l){n + l) \ J ) 

where 

{z,x} = (z'z<® - (3/2)z" 2 ) z'- 2 (4.8b) 

is the Schwarzian derivative, and where z — dz/dx, we can put Eq. (|2.3|) into a 
form in which the terms of orders n — 1 and n—2 do not appear. After the renaming 
of variables and coefficients of the new equation using the same notation as for the 
original equation 1|2.3|) . the transformed equation takes the form 



(n) + a 3 j/"- 3) + • • • + a n y = 0. (4.9) 
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In fact, due to the prominence in size of invariants of differential equations and 
the rapid rate at which this size increases with the number of nonzero coefficients 
in the equation, it is customary to use the canonical form (|4.9I) of (|2.3p for the 
investigation of invariants [6] [3] . 

A determination of invariants of linear ODEs by means of the infinitesimal gen- 
erator of the form (|4.6p was made in [14] for equations of order up to five in the 
various canonical forms (|2.3[1 . (|4.7|) . and (|4.9|) . but only for low orders of prolon- 
gation of the operator X a not exceeding three. We now undertake the explicit 
determination of a generating system of invariants of all orders by an application 
of formulas (|4.3[) - (|4.5[) . For this purpose we shall adopt the canonical form (|4. 9|) for 
which the invariants have a much simpler expression, and we restrict our attention 
to the case of equations of order n = 5 which have the form 

y i5) + a 3 y" + a iV ' + a 5 y = 0. (4.10) 

For every n > 3, and for any order p of prolongation of the operator X°, the number 
F of absolute invariants of Eq. (|4.9|) can be shown [14] to be given by 

r =n + 4-p(n-2), (4.11) 

and this indicates that for a fixed value of n, as the order of prolongation of X° 
increases by one unit, the number of absolute invariants increases by n — 2, which is 
precisely the number of coefficients in the equation. This means that a generating 
system of absolute invariants should contain n— 2 elements, which by virtue of (|4.11[) 
can be taken to be functionally independent, provided that none of the coefficients 
cij vanishes. For (|4.10j) . X° depends on three arbitrary constants k\, k 2 and k 3 and 
has the explicit form 



(4.12) 



X° = [ki + x(k 2 + k 3 x)] d x - 3a 3 (k 2 + 2k 3 x) d a3 

- 2 [-3a 3 fc 3 + 2a 4 (fc 2 + 2k 3 x)] d ai + [-4a 4 fc 3 - 5a 5 (fc 2 + 2k 3 x)] <9 Q5 

and its first prolongation has four invariants given by 

h jl^L (4 , 3a) 

= (Z21+4)! (4 . 13b) 

216 a 3 

5a| + 9aja' 5 - 3a 4 a 3 (5a 5 + 2a' 4 ) + 3a| a 3 

/3 94 (4 ' 13d) 

The fundamental relative invariant of (|4.9p is readily seen to be given by So = a 3 (x), 
and Eq. (|4.13|) shows that the four functions 

Ro = 3a 5 a 3 — a 4 

Si = — a 4 + a 3 

5*2 = 6a 3 a 4 — a 4 — 6a 4 a 3 

S3 = 5a 4 + 9a 3 a' 5 — 3a 4 a 3 (5a 5 + 2a' 4 ) + 3a 4 a 3 
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are relative invariants of respective index 8,4,8, and 12. In terms of these relative 
invariants, we have 

Ia = ~oi' h = cT' 12 = of' ^ 3 = 73- (4-14) 
o o o j 

Moreover, each of the four functions Ro,Si,S2, and 53 is of order at most one, and 
together they form a functionally independent set. Consequently, on in accordance 
with (|4.5jl . B = {Io,Ii,l2,h} forms a functionally independent generating system 
of absolute invariants which may be called a basis of absolute invariants of (|4.9|) . 

It also follows that a fundamental set of absolute invariants of order p > 1 is 
given by the functions 

h, Xk(Si), Xk(S 2 ), Xk(S 3 ), iork = 0,...,p-l (4.15) 

and their number is 3p + 1, which is in accordance with equation (|4.11[) . It follows 
from Eq. (|4.5[) that to find explicit expressions for the absolute invariant Xk(Sj), we 
only need to compute the corresponding relative invariant tpk(Sj) and to determine 
its index. For Si, the corresponding higher order relative invariants up to the order 
three are 

(p(Si) = 4(a 4 - a 3 )a 3 + 3a 3 (-a 4 + Og) 
l P2 (S 1 ) = 32(a 4 - a 3 )a 3 2 - 3a 3 (9a 4 a 3 + (4a 4 - 13a 3 )a 3 ) 
+ 9al(a'l-4 3y ) 

and their indices are clearly 8 and 12, respectively. The size of these two invari- 
ants are much more smaller compared to those corresponding to S2 and S3 , simply 
because Si has a much smaller size compared to S2 and S3 . For instance the expres- 
sion for <fik{S3) is about five times larger in size than that for ipk(Si), for A; = 1,2. 
It turns out that the expression for the invariants computed explicitly by means 
of the infinitesimal generator X° are considerably much smaller in size than those 
computed from the indefinite sequence (|4.5[> . Indeed, a direct computation shows 
that a fundamental set of differential invariants of order two is given by 

(7a\ - 14a 4 a 3 + 6a 3 a 3 ) 3 

la — o 

216 a% 

_ 4a| + 24a|a 3 + 9a|a 4 ' - 9a 4 a 3 (3< + a 3 ') 

15 = ' n 

9d3 

_ (-18a| - 18a|a 3 + 18a 3 a' 5 ' - 6a^al(lla' 5 + 2a'() + a|a 3 (55a 5 + 40< + 6a 3 ')) 3 



5832 a 



and that for differential invariants of order three is given by 



I- 
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-14a| + 42a^a 3 - 36a 4 a 3 a 3 ' + 9a§a 3 3) 



9 a| 



(-2ai - 12a|a 3 + 3ala 3 (5a 4 + 3a 3 ') + - 2a 4 a§(5a 4 ' + a 3 3 } )) 3 



8al 6 



h = (5 9il + S 9 , 2 ) 3 /(5832 a f) 
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where 

5 9 .1 = 35a| + 45a|a 3 - 10ata 3 (lla 5 + lla' A + 3a 3 ') + 18a^4 3) 

59.2 = -12a 4 a|(9ag + af ) + 6^(334 + lla" + 4 3) ) 3 ). 

As already mentioned these absolute invariants of (|4.10[) are given in [14] , but only 
for the second prolongation of the operator X , and not in a form in which the index 
can be readily read off. They indeed appear to be smaller in size than those derived 
from the recurrence relations (|4. 5[) , which indicates that a further simplification of 
the expression of the function if might be possible. However, such a simplification 
might lead to more complicated recurrence equations for the indefinite sequence of 
invariants. It's however naturally much easier to generate higher order invariants 
using the indefinite sequence (|4.5|) . rather than finding first the right prolongation 
of X° and then solving the corresponding system of equations to find the invariants. 

Similarly, using Eqs. (|4.14[) , (|4. 13[) . and (|4.5[) . we can also determine, by invoking 
Corollary dj a fundamental set of relative invariants of all orders. For instance, 
fundamental relative invariants of order up to two of (14. 10|) are given by 

<Pk{Sj,S ), j = 1,2,3, fe = 0,l (4.16a) 
WkCRo.Sb), Vk(S ,Ro), fc = 0,1,2. (4.16b) 

Incidently, Halphen [5] constructed by a different method an indefinite sequence 
of relative invariants similar to that given in (|4.16|) , but for linear ODEs of order 
4 in the canonical form (|2.3p . In particular the approach of [5] is not based on 
infinitesimal nor Lie groups methods and it aimed only at deriving and indefinite 
sequence of relative invariants from known ones. Consequently, the sequence of rel- 
ative invariants obtained in the said paper for equations of order 4 is not associated 
with the determination of a fundamental set of absolute invariants. 
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